Single-shot Sub-Rayleigh Imaging with Sparse Detection 
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For conventional imaging, the imaging resolution limit is given by the Rayleigh criterion. Ex- 
ploiting the prior knowledge of imaging object's sparsity and fixed optical system, imaging beyond 
the conventional Rayleigh limit, which is backed up by numerical simulation and experiments, is 
achieved by illuminating the object with single-shot thermal light and detecting the object's informa- 
tion at the imaging plane with some sparse-array single-pixel detectors. The quality of sub-Rayleigh 
imaging with sparse detection is also shown to be related to the effective number of single-pixel 
detectors and the detection signal-to-noise ratio at the imaging plane. 
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For standard conventional imaging of directly record- 
ing the object's intensity distribution with a charge cou- 
pled device (CCD) camera, both the imaging system's 
Rayleigh limit and the camera's pixel-resolution restrict 
the optical system's imaging resolution [l|. For example, 
the imaging resolution is mainly determined by the op- 
tical system's Rayleigh limit in remote sensing because 
the numerical aperture {N.A.) of imaging lens is usually 
small relative to the detection distance. While in long- 
wavelength radiation band such as infrared and terahertz 
imaging, because the camera with large planar arrays 
is very hard to manufacture, the imaging resolution is 
mainly limited by the camera's pixel-resolution. 

Compared with overcoming the camera's pixel- 
resolution to imaging resolution, many methods are in- 
vented to overcome the imaging system's Rayleigh limit 
at present 0-[l3]- Exploiting the evanescent components 
at the object's immediate proximity, sub-Rayleigh imag- 
ing can be achieved, but this method is only applied 
in the near-field range 0-01 • Several microscopy tech- 
niques based on fluorescence are also introduced to im- 
prove imaging resolution. However, it requires scanning 
or repetitive experiments, thus which limits real-time ap- 
plications [^-Q. In addition, using additional a priori 
information of optical system, the imaging resolution be- 
yond Rayleigh diffraction limit can be obtained. How- 
ever, the improvement degree is limited in practice be- 
cause of the influence of detection noise [l0l - [l4{ . 

Recently, the image's sparsity has been taken as popu- 
lar a priori, which is a quite general assumption because 
most of the natural objects are sparse in a known basis 
(or under a suitable basis transform), an image can be 
stably reconstructed by sparse reconstruction technique 
even if the measurement number is less than Nyquist rate 
and this technique, to some extent, is robust to noise [15- 
[igj . Moreover this technique has already been applied 
to de- noising [ll. Ell , super-resolution imaging [20, |2l| , 



remote sensing [22, |23| , compressive imaging |24|-|28| , and 
magnetic resonance imaging with success [2l|. For con- 
ventional imaging, the optical system is usually fixed 
or known, so the optical system's point spread function 
(PSF) can be also taken as a priori. Furthermore, based 
on the object's sparse priori property, measurements be- 
low Nyquist rate are required to exactly restore the ob- 
ject. Therefore, exploiting the object's sparsity assump- 
tion in a known basis and the prior knowledge of PSF, 
single-shot sub-Rayleigh imaging is possible even if some 
sparse-array single-pixel detectors are used to record the 
object's information at the imaging plane. 
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FIG. 1: Schematic of thermal-light single-shot imaging with 
sparse detection. 

Fig. 1 presents experimental demonstration schematic 
of thermal-light single-shot sub-Rayleigh imaging with 
sparse detection. The uniform light emitting from a halo- 
gen lamp is filtered by an optical filter (with the center 
wavelength A=650 nm and the bandwidth AA=10 nm) 
and then collimated by a lens with the focal length /q. 
The object is illuminated by the collimation light and 
then imaged onto the imaging plane xt by a standard 
conventional imaging setup. Different from recording the 
object's transmission information by a CCD camera at 
the imaging plane Xt, the detection system Dt we pro- 
pose is characterized by some sparse-array single-pixel 
detectors. 

The distances zi , Z2 and the focal length of the lens / 
obey Gaussian thin-lens equation: + — j- Based 
on Rayleigh criterion [l[ , the resolution limit Aa;^ of con- 
ventional imaging shown in Fig. 1 is determined by the 
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wavelength A and the N.A. of the lens /, namely 

Xzi 



Axs = 0.61 



A 

N.A. 



1.22- 



(1) 



where L is the effective transmission aperture of the 
imaging lens / and A'^.^. is approximate to j^. 

When the light is fully spatially incoherent and uni- 
form, the intensity at the imaging plane xt is the con- 
volution of the intensity at the object plane with the 
incoherent point spread function (PSF) [l3| 



Itix) = lobjix) ® h{x) + Inoiseix). 



(2) 



where h{x) is the optical system's PSF, Iobj{x) is the 
object's intensity distribution, Inoise{x) is the intensity 
distribution of the noise and It (x) is the intensity distri- 
bution at the imaging plane xt. For the optical system 
shown in Fig. 1, the PSF is 



h{x) = sinc^[-- — x] 

Xzi 



i(7ra;) 



(3) 



where sin c(x) 

From Eqs. (2) and (3), compared with the original 
object, it is a blurred image with low spatial resolution 
at the imaging plane Xt when the transmission aperture 
of the imaging lens / is small. In order to restore a 
high-resolution image, the relation shown in Eq. (2) 
usually yields F{It{x)} = T{Iobj{x)}H + I'{Iobj{x)} in 
the spatial-frequency domain, where J- denotes Fourier 
transform and H is the optical transfer function. There- 
fore, the problem is that we wish to reconstruct the ob- 
ject but the measured image is smeared by a low-pass 
filter. Although many iterative methods were used to re- 
store the object, the reconstruction quality strongly de- 
pends on the signal-to-noise ratio (SNR) in the measured 
data Furthermore, for the optical system shown 

in Fig. 1, because the detection system Dt is charac- 
terized by some sparse-array single-pixel detectors and 
many low-frequency information will be also missed, it 
becomes much more difficult to restore the object by tra- 
ditional iterative methods, compared with recording all 
the smeared image's low-frequency information by a CCD 
camera. 

Different from the method described above, we try to 
directly restore the object borne on the optical system 
shown in Fig. 1, exploiting the sparsity assumption of the 
object in a known basis and the prior knowledge of the 
fixed optical system. According to sparse reconstruction 
theory, even if the measurement process is noiseless, there 
are an infinite number of images, which— after being con- 
voluted by the PSF— will result in the smeared image, our 
convex optimization program is how to find the sparsest 
one. The sparse reconstruction technique has mathemat- 
ically demonstrated that if the object is sparse enough, 
then any sparsity-based reconstruction method is bound 
to find the sparsest solution [l5l - [l9| . Here, we have em- 
ployed the gradient projection for sparse reconstruction 



algorithm '30|, the object T can be reconstructed by solv- 
ing the following convex optimization program: 



/ 1 2 

T = \t'\; which minimizes : - It{xi) — h{x) ® \t' {x)\ 



+T 



^{\t'{x)\'} ^,V, = l.-.Af.(4) 



where r is a nonnegative parameter which is used to ad- 
just the sparse weight in the image restoration process, 
M is the number of single-pixel detectors at the imaging 
plane Xt and ^ denotes the transform operator to the 
sparse basis. WVW^ and \\V\\-y^ represent the Euclidean 
norm and the ^i-norm of V ^ respectively. 
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FIG. 2: The simulated demonstration of noiseless sub- 
Rayleigh imaging with sparse detection, (a) The original 
object; (bi) the intensity distribution recorded by M=4096 
sparse-array single-pixel detectors at the imaging plane xt and 
the distance between two singe-pixel detectors in lateral di- 
rection Aa;=6.45 /im; (b2) M=1024 and Aa;=12.90 ^im; (bs) 
M=256 and Aa;=25.80 ^im; (b4) M=64 and Aa;=51.60 ^im; 
(bs) M=36 and Aa;=64.50 ura; (be) M=72 and Aa;=64.50 
/im; (ci-ce) are corresponding sparse reconstruction results 
with respect to (bi-bs). 

In order to verify the idea, Fig. 2 and Fig. 3 have 
given the simulated and experimental demonstration of 
sub-Rayleigh imaging with sparse detection, using the 
optical system depicted in Fig. 1. The original object 
(64x64 pixels, and the pixel size is 6.45 /imx6.45 /im), 
as shown in Fig. 2(a) and Fig. 3(a), is a double-slit 
(slit width a=30 /xm, silt height h=240 /im and center- 
to-center separation d=60 /xm) and it comprises of 370 
nonzero values in real-space domain. The parameters 
listed in Fig. 1 are set as follows: zi=Z2=800 mm, the 
focal length of the lens /=400 mm and its effective trans- 
mission aperture L=8.0 mm. In addition, the pixel size 
of the single-pixel detector is 6.45 /imx6.45 /im. There- 
fore, based on Eq. (1), the imaging system's resolution 
limit is Axs ~80 /im and the object's image at the imag- 
ing plane Xt can not be resolved for conventional imag- 
ing [Fig. 2(bi), Fig. 3(bi)]. However, both the simu- 
lated and experimental results illustrated in Fig. 2 and 
Fig. 3 clearly demonstrate that sub-Rayleigh imaging 
with sparse detection can be achieved by exploiting the 
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FIG. 3: The experimental demonstration of sub-Rayleigh 
imaging with sparse detection, (a) The original object im- 
aged by a conventional optical imaging setup with large N.A.; 
(bi-be) and (ci-ce) are the intensity distributions record by 
the detectors at the imaging plane xt and their corresponding 
sparse reconstruction results, the same as Fig. 2. 
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FIG. 4: Sub-Rayleigh imaging with sparse detection in dif- 
ferent detection SNR and (M, Ax). (a) The intensity 
distributions recorded by M=1024 (Aa;=12.90 /im) sparse- 
array single-pixel detectors and the detection SNR=2 dB 
at the imaging plane xt; (b) M=1024 (Aa::=12.90 /im) and 
SNR=5 dB; (c) M=1024 (Ax=12.90 fim) and SNR=10 dB; 
(d) M=1024 (A2:=12.90 fim) and SNR=15 dB; (e) A/=120 
(Aa:=38.70 /xm) and SNR=15 dB; (f) M=72 (Aa;=64.50 fim) 
and SNR=15 dB; (ai-fi) are corresponding sparse restoration 
results with respect to (a-f). 

prior knowledge of the object's sparsity and the optical 
system's PSF. When the number of singe-pixel detectors 
(namely the total effective measurement number M) at 
the imaging plane xt is M=4096, 1024, 256, 64, 36, and 
72, Fig. 2(bi-b6) and Fig. 3(bi-b6) present the intensity 
distributions recorded by the detectors, and their corre- 
sponding sparse reconstruction results borne on Eq. (4) 



are depicted in Fig. 2(ci-C6) and Fig. 3(ci-C6), respec- 
tively. As the distance between two singe-pixel detectors 
is increased, the restoration quality will be reduced be- 
cause of the decrease of the sampling number [Fig. 2(ci- 
C5) and Fig. 3(01-05)]. However, if the sampling number 
at the imaging plane Xt is increased, sub-Rayleigh imag- 
ing with sparse detection can be stably restored, even if 
the distance between two singe-pixel detectors (Aa;=64.5 
fim) is large than the double-slit's center-to-center sepa- 
ration [Fig. 2(05-05) and Fig. 3(05-05)]. In addition, as 
shown in Fig. 2(05) and Fig. 3(05), only M=72 measure- 
ments are used to restore the object's image, which is far 
fewer than suggested by Shannon's sampling theorem. 

As shown in Fig. 4, we also perform the dependance 
of sub-Rayleigh imaging with sparse detection on the 
detection SNR at the imaging plane Xt- Fig. 4(a-f) 
present the intensity distributions recorded by the de- 
tectors at the imaging plane xt in different SNR. Based 
on the sparse reconstruction technique described by Eq. 
(4), the corresponding restoration results are displayed 
in Fig. 4(ai-fi). From Fig. 4(ai-di), even if the de- 
tection SNR is 2 dB, we can approximately reconstruct 
sub-Rayleigh imaging using M=1024 measurements and 
the reconstruction quality will be improved as the in- 
crease of the detection SNR. While in the case of the 
detection SNR=15 dB, high-quality sub-Rayleigh imag- 
ing with sparse detection, as shown in Fig. 4(di-fi), can 
be still restored by using the measurements far below 
Nyquist rate. 

In conclusion, we have realized single-shot sub- 
Rayleigh imaging with sparse detection, exploiting the 
imaging object's sparsity and the prior knowledge of fixed 
optical system. Both the simulated and experimental 
results have demonstrated that using the measurement 
below Nyquist rate, sub-Rayleigh imaging with sparse 
detection can break through the limitation of both the 
optical system's Rayleigh limit and the camera's pixel- 
resolution to imaging resolution. Furthermore, the ex- 
periment has also shown that sub-Rayleigh imaging with 
sparse detection can still be approximately reconstructed 
in the case of measured SNR=2dB. This technique is very 
useful to microscopy of living cells or bacteria, and imag- 
ing of atoms captured by ion trap, etc, where the images 
are enough sparse. 
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